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We experimentally demonstrate how magnetic depolarization of velocity-selective optical pump- 
ing can be used to single out the collisional cusp kernel best describing spin- and velocity relaxing 
collisions between potassium atoms and low pressure helium. The range of pressures and transverse 
fields used simulate the novel optical pumping regime pertinent to sodium guidestars employed in 
adaptive optics. We measure the precession of spin-velocity modes under the application of trans- 
verse magnetic fields, simulating the natural configuration of mesospheric sodium optical pumping 
in the geomagnetic field. We also provide a full theoretical account of the experimental data using 
the recently developed cusp kernels, which realistically quantify velocity damping collisions in this 
novel optical pumping regime. A single cusp kernel with a sharpness s — 13 ± 2 provides a global 
fit to the K-He data. 
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Wc demonstrate a new way to measure the transition 
matrix W(v, v') for velocity-changing collisions between 
ground-state alkali metal atoms and arbitrary buffer gas 
molecules in the classical limit. This matrix contains all 
the information required for computing the contribution 
of collisions to the evolution of non-equilibrium velocity 
distributions in the vapor, which is given in general by 



7vdJ>(«,*) + 7vd / dv'W(v,v')p(v',t) (1) 



dp(v, t) 
dt 



where p(v,t) is the time-dependent velocity distribution 
and 7 v d the rate of velocity-changing collisions. Precise 
measurements of this matrix have long been of inter- 
est as a way to place empirical constraints on the in- 
teratomic potential for alkali metal - buffer gas pairs 
[l[. Two new applications demand an increasingly pre- 
cise knowledge of W(v,v'). First of all, the ease with 
which non-equilibrium velocity distributions can be cre- 
ated through velocity selective optical pumping makes 
alkali-metal vapor cells an ideal candidate for investigat- 
ing new theorems of non-equilibrium statistical mechan- 
ics ,v') is the transition probability that arises 
in Markov derivations of these theorems, and quantita- 
tive knowledge of its behavior for real systems would 
make such experiments easier to design and analyze. 
Secondly, laser guide stars employed in the adaptive 
optics systems of modern terrestrial observatories pro- 
duce non-equilibrium velocity distributions of optically 
pumped mesospheric sodium atoms. The numerical mod- 
els used for optimizing the laser parameters for maximum 
backscattcr require W(v, v') as input, but the dearth 
of empirical measurements for sodium-air collisions has 



forced current models to use a highly simplified model, 
where W(v,v')is the Boltzmann distribution in v inde- 
pendent of v' Q . 



'Present address: Department of Physics, Columbia University, 
538 West 120th Street, New York, NY 10027-5255, USA. 



I. THEORY 

The theoretical framework we employed to design the 
experiment and analyze our data is laid out in detail in 
U . Here we summarize the points that are most directly 
relevant to the current work. 

We model the vapor cell as an ensemble of single 
atoms whose velocities are perturbed by collisions with 
molecules of an inert buffer gas and with the cell walls. 
Neglecting radiation-pressure effects relevant to laser 
cooling experiments, the equilibrium velocity distribution 
of atoms in the buffer gas is the Maxwell-Boltzmann dis- 
tribution. It will be convenient to measure atomic veloci- 
ties in units of the most probable speed Vp = y / 2ksT/M 
where T is the temperature of the cell and M is the 
mass of the atom. We will denote the projection v of 
the atomic velocity onto the laser beam direction by 
x = v/vd- Although the overall velocity distribution 
remains constant, the incident laser beam introduces 
correlations between the spins and the velocities: due 
to the Dopplcr shift, a laser at frequency ui excites a 
resonance at frequency too only in atoms with velocity 
.To = (lu — LU(j)/kvD along the beam. This is the origin 
of 'velocity-selective optical pumping' (VSOP), whereby 
the distribution of atomic spin states becomes velocity- 
dependent. We assume the laser is weak enough that the 
small fraction of atoms in the excited states at any given 
time can be ignored, so the overall effect of the laser is to 
modify the distribution of ground-state sublevels at the 
resonant velocities. The distribution among energy sub- 
levels \/j.) for atoms with velocities between x and x + dx 
is described by the (2/+l)(2S + l) x (2/+l)(25 + l) den- 
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FIG. 1: (Color online) (a) The monochromatic, circularly po- 
larized pump laser of frequency uj selectively spin-polarizes 
atoms in the velocity group around that is Doppler-shifted 
into resonance. The counter-propagating probe beam of the 
same laser probes atoms at the velocity group around — v u . 
Pumped atoms collide with buffer gas atoms and diffuse in 
velocity space from the velocity v' before to the velocity v 
after the collision. The velocity v follows the probability dis- 
tribution W(v,v'), shown in (b) for v' = w w . By scanning 
the laser frequency to the probe beam measures W(— v u , v u ). 
A transverse magnetic field causes the spins to precess dur- 
ing this process, acting as a "clock" during the velocity-space 
diffusion in between pumping and probing, (c) The overall 
effect of all the collisions is to produce a steady-state dis- 
tribution of average spin as a function of velocity given by 
W(v,Vu,) for two different values of the transverse magnetic 
field B± . The effect of spin precession is seen in the sign re- 
versal of W{— Vu, v u ) at B± = 0.5 G relative to W(— v^jV^) 
at B± = G. 



sity matrix dp(x)dx = Xtiv{ x )\^) { v \ dx, where S and 
/ are the electron and nuclear spins, respectively, and 
and \v) are eigenstates of the ground-state Hamiltonian. 
The latter includes hyperfinc splitting and Zccman cou- 
pling to an external magnetic field and the corresponding 
eigenvalues are and E v . 

In a VSOP situation with no velocity-changing colli- 
sions, some of these modes will have distributions X^i x ) 
sharply peaked around the resonant velocities xq. As 
illustrated in Fig. [TJ adding an inert buffer gas al- 
lows atoms pumped at the resonant velocities to jump 
to other, non-resonant velocities through collisions with 



the buffer gas molecules that change the atomic velocity 
while preserving spin polarization. Fig. [T] also displays 
the effect of a transverse magnetic field, namely the spin 
will undergo Larmor precession while the atom diffuses 
from the pump resonant velocity to the probe resonant 
velocity, thus adding a new tunable timcscalc to the prob- 
lem. 

The velocity-changing collisions are characterized by 
the 'collision kernel' W(x,x'), which gives the probabil- 
ity that an atom with initial velocity x' will have veloc- 
ity x after a single collision. Since our atoms remain in 
the spherically symmetric S state of orbital angular mo- 
mentum for the overwhelming majority of the time, the 
collision kernel does not bear any dependence on the in- 
ternal atomic state. The overall effect of collisions is to 
give Xfj.v( x ) a broad pedestal around the sharp resonant 
spike, and the shape of this pedestal contains valuable 
information about the collision kernel. To extract this 
information, we need to first write the steady-state solu- 
tions to the evolution equation for the X^'s in terms of 
W(x,x'), and then find a good ansatz form for W(x,x') 
that can reproduce experimental data through appropri- 
ate adjustment of a few free parameters. 

The evolution of the velocity distributions Xiiv 

(x) of 

the spin modes is described by 

^^dt' ^ = ~~ J x ')X^u{x' , t)dx' + P^y (x). (2) 

where the damping kernel is 

K fiv (x,x') = [7 w +i^+7vd]£(x--:r')-7 vd VF(2;,x'). (3) 

The rate of loss of polarized atoms due to diffusion or 
free flight to the walls is 7 W . The Bohr frequencies 
fl^ v = (E^ — E u )/H depend on the applied magnetic 
field. For laboratory experiments with alkali-metal atoms 
in low buffer-gas pressures, the direct damping of spin 
modes by gas-phase collisions can be neglected, since 
practical buffer gases like N 2 or the noble gases He, Ne, 
Ar, etc. have such small spin- flip rates compared to j w . 
The main effect of the buffer gas on spin relaxation is to 
lengthen the time needed for a polarized atom to diffuse 
to the wall, that is, to decrease 7 W . The rate of velocity 
relaxation due to the buffer gas is parametrized by the 
rate 7 V( j at which atoms with initial velocity x' are trans- 
ferred to a range of final velocities x. The distribution of 
final velocities is given by the collision kernel W{x, x'), 
which is the function we are trying to measure. 

The laser-atom interaction that creates the non- 
equilibrium distribution of optically pumped atoms 
enters the model through the source term P^ v = 
P^vix^x^'v')-! which depends on the values x^'v'ix) of 
all the spin mode amplitudes at velocity x. For a true 
monochromatic light source exciting electrons to states 
with infinite lifetimes, P^ u would be a superposition of 
delta functions. In fact, the finite excited-state lifetime, 
finite laser linewidth and other experimental issues to be 
discussed later combine to produce a Lorentzian source 
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term with a particular width, which can be empirically 
determined from a measurement without buffer gas. 

The noncquilibrium steady-state velocity distribution 
of the spin modes, found by setting dxuv/dt = in Eq. 
([5]), is formally given by 



X^(x)= K-^{x,x')P^{x')dx', 



(4) 



where the inverse damping kernel is defined by 
/ K~^(x,x')K^ v (x',x")dx' = S(x - x"). Note that 
depends on the x^'s, so this is not a closed- form solu- 
tion, but we can use it to obtain a power series solution 
in the pump laser intensity, as shown explicitly in Eqs. 
(134) through (136) of 0|. Our numerical calculations 
use the first-order solution, given by Eqs. (137) through 
(140) of 0| . 

We now need a physically reasonable and computa- 
tionally tractable Ansatz for W(x, x') in order to extract 
information about the collisions from experimental sig- 
nals. McGuyer et al. have shown that 'cusp kernels' 
C s (x, x') meet these requirements @. The cusp kernel 
function has a single free parameter s called the 'sharp- 
ness,' which we will employ as a fit parameter to match 
our model predictions to our experimental data. Setting 
W(x, x') = C s (x,x'), the inverse damping kernel K~l 
takes on the particularly simple form 



K-l{x,x') 



1 



6(x - x') + 



7w + i^nv + 7vd 

7vd 



7 W + iil 



I"' 



where 



\1V 



7w ~l~ iQp,v 
7 W + ifi M „ + 7 v d 



s. 



(5) 



(6) 



Using this result and Eq. (|4]) , we can obtain the steady- 
state Xp.v ' s m terms of the atomic parameters, 7 W , 7 v d 
and s. The rates 7 W and 7 V( j can be determined experi- 
mentally, so the solution has s as its only free parameter. 
Once we calculate the probe beam absorption spectrum 
from the Xnv' s > we can easily find the sharpness s that 
gives a least-squares best fit to the data. 

The experiment we performed is similar to the exper- 
iment of Aminoff et al. @, the main difference being 
the addition of a transverse magnetic field introducing 
the magnetic depolarization discussed previously. The 
authors in Q noted that the traditionally used Keilson- 
Storer kernel was a poor approximation of the actual ker- 
nel, since it could not fit the data over a broad range of 
pressures and the kernel's parameter was thus not unique 
but pressure-dependent. In contrast, we will demonstrate 
that a cusp kernel with a single sharpness can provide a 
global fit to the data. The additional, experimentally 
controllable parameter of the transverse magnetic field, 
incorporated in our formalism through the frequencies 
f2^„, provides a systematic check of our understanding 
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FIG. 2: (Color Online) Counter-propagating pump-probe 
lasers create and probe spin-velocity distributions of potas- 
sium metal in helium buffer gas. The atoms are contained in 
sealed glass cells of 1/2 in. diameter and 5/4 in. length, re- 
siding in an oven air-heated at 69 ±0.2 °C to maintain a high 
enough potassium vapor pressure of natural isotopic abun- 
dance, 93.26% 39 K and 6.73% 41 K. The cell and oven reside in 
the middle of three pairs of Helmholtz coils along three mutu- 
ally orthogonal axes providing a homogeneous, constant and 
tunable magnetic field. The intensities of the two counter- 
propagating beams were independently controlled with two 
neutral density filters. The beams are not perfectly antipar- 
allel due to the pump mirror blocking the probe detection 
photodiode, but the angle between them was less than 5 mrad. 



of the kernel, as it introduces an additional degree of 
freedom in the spin-velocity time evolution, not present 
in the previous experiments using longitudinal magnetic 
fields only. 



II. EXPERIMENT AND DATA ANALYSIS 

In. Fig. [2] we depict the schematic of our experiment. 
We inserted three cylindrical sealed glass cells containing 
potassium and helium buffer gas at three different pres- 
sures into the path of two counter-propagating circularly 
polarized beams from a Toptica DL Pro tunable diode 
laser. The pump beam is responsible of VSOP, while 
we monitor the transmission of the weaker probe beam. 
To amplify the signal, we modulated the circular polar- 
ization of the probe beam with a photo-elastic modula- 
tor (PEM), cycling the polarization from right-circular 
to left-circular and back at a frequency of 42 kHz. Since 
both polarizations should have the same transmission ef- 
ficiency with no optical pumping, we were able to extract 
and amplify the correction to the lineshape due to opti- 
cal pumping using a lock-in amplifier phase-locked to the 
PEM, while slowly scanning the laser frequency over the 
potassium's Doppler profile. 

Since the pump beam only pumps atoms having veloc- 
ities x^c — (to — Wjkp)/ kvo j where w^p is the frequency 
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p He = 32 mtorr 
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FIG. 3: (Color) Potassium-Helium data. Shown is the measured (red) and calculated (blue) absorption rate of the probe beam 
for three different magnetic fields (zero field, 0.2 G and 0.5 G, transverse to the laser propagation direction) at pressures of 32, 
65, and 120 mTorr. The theoretical calculation was fit to the data by a least squares fit. A single cusp kernel with sharpness 
s = 1 3 provides a good fit to all data. 



of the transition from ground state sublevel /i to excited 
state sublevel i>, and the probe beam only interacts with 
atoms with opposite velocities — the spectrum of 
probe beam transmission will display sharp saturated- 
absorption spikes at uj^ as well as the cross-over lines 
(uj^p + uj^,-/)/2, where the pump and probe beams can 
simultaneously interact with the same atoms. To ensure 
that the probe beam does not modify the atomic den- 
sity matrix, we measured the height of the tallest one 
of these spikes at various probe intensities and a con- 
stant pump intensity. At low probe intensities, the height 
scales linearly with intensity, while it becomes quadratic 
as the probe intensity is increased. By chosing a ratio 
of probe to pump intensity of about 10%, we ensured we 
are well within the linear regime. Similarly, we can ne- 
glect the second and higher-order terms in pump beam 
intensity from the series solution of Eq. @ by examining 
the resonance's scaling with pump intensity and staying 
within the linear regime. As indicated in Fig. [2l the 
filter (NDF1) for adjusting the pump intensity was posi- 
tioned before the beamsplitter, so that we could adjust 
the pump intensity at constant probe/pump intensity ra- 
tio. In the first-order regime, the transmitted intensity 
should scale linearly with each beam intensity individu- 
ally, and so their simultaneous variation should produce 
a quadratic shape. We set NDF1 so that we are well 
within this quadratic regime. 

In Fig. |3] we show the data acquired from our three 



cells at three different transverse magnetic fields for each 
cell, along with the theoretical fits. To compare the mea- 
sured lineshapes to the cusp kernel predictions we used 
Eqs. (143) through (147) of Q to obtain the first-order 
lineshape corrections from the spin mode amplitudes Xp.v- 
The required rates 7 v d and 7 W were obtained from the ex- 
pressions 



and 



7w 



7vd 



a 

vd 



2zr 



a 

8D 
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(7) 



(8) 



given in Eqs. (Al) and (79) of Q respectively. The most 
probable speed vd and the sharpness parameter s have 
already been defined. The laser beam and cell radii are 
a = 0.06 cm and b — 0.56 cm, respectively, and D is 
the diffusion coefficient of potassium atoms in helium, 
inversely proportional to the helium pressure. We used 
the value Do = 0.45 cm 2 /s for 1 atm of helium. Note 
that Eq. {?]) is valid if the collision kernel is adequately 
modeled by a single cusp kernel. 

To correct for potassium's optical thickness, leading to 
pump and probe intensity variation along the length of 
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the cell, we divide the data by the factor 

/M = -^-(l-e-^) (9) 

where <t(u) is the absorption cross-section of the un- 
polarizcd potassium atoms, n is the number density of 
potassium atoms (a known function of temperature) , and 
I = 2.1 cm is the length of the cell. 

We also had to account for the fact that our pressure 
gauge had a large uncertainty at the low pressures we 
were using. To get a good fit, we had to adjust the pres- 
sure values input to the model within the range of un- 
certainty, and ultimately used the values 32, 65, and 120 
mTorr instead of the nominal 25, 50, and 100. Finally, 
we added one more free parameter to match the absolute 
scale of the y-axis in Fig. [3]between model and data. We 
used the same scale factor for all nine data traces in Fig. 
[3l The simultaneous least-squares fit of the model to the 
nine data traces is superimposed on the data in Fig. [31 
resulting in a sharpness of s = 13 ± 2. 



III. CONCLUSIONS 



In conclusion, we have demonstrated a simple way 
to accurately determine the sharpness of cusp kernels 
modeling velocity-selective optical pumping. In par- 
ticular, we have shown that a single cusp kernel is 
sufficient to reproduce lineshape measurements far from 
the Maxwell-Boltzmann equilibrium, across a variety of 
pressures and magnetic fields. The employed method of 
magnetic depolarization allows a rapid, experimentally 
straightforward and accurate measurement of collision 
kernels useful for constraining interatomic potentials, 
testing non-equilibrium statistical mechanics theorems, 
or optimizing laser guidestars systems. 
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